function MU1 = MU1(c,gamma,cf)

MU1=cf*c^(-gamma);


end